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DOPPLER FREQUENCY IN INTERPLANETARY RADAR 
AND GENERAL RELATIVITY 


by 

G. C. McVittie 

University of Illinois Observatory 


INTRODUCTION 


The object of this investigation is to present the theory of the change of frequency of signals in 
interplanetary radar experiments when general relativity is employed as the theory of gravitation. The 
basic ideas will be found in the author’s General Relativity and Cosmology (Reference 1 ). No claim is 
made that the contents of the first three sections of this paper are entirely hew and that similar results 
have not been obtained by others by alternative methods. The last two sections contain the author’s 
ideas on the interpretation of the theoretical formulas. The associated problem of radar time delays is 
worked out in Reference 2. 

It is well known that in Einstein’s theory of gravitation an exact solution of the gravitational field 
of a set of discrete bodies is possible only when one of the bodies is of finite mass whereas the rest are 
of infinitesimally small mass. This is in contrast to Newton’s theory of gravitation in which an exact 
solution for the problem of two massive bodies in motion is possible, complications arising only when 
three or more bodies are in question. In Einstein’s theory the field of the single massive body is given 
by the Schwarzschild space-time. The massive body may be taken to be the Sun and the planets, their 
satellites, space probes, etc., must then be regarded as test particles whose own gravitational fields are 
to be regarded as negligibly small. 


The simplest form of the metric of the Schwarzschild space-time is 

, / 2 m\ , dr 2 /( 1 - 2m/r) + r 2 dCl 2 

ds L = ( 1 - — )dr 

V ' / c 2 


( 1 ) 


dn 2 = dd 2 + sin 2 e dtp 2 , 

in which r = 0 is the center of the (spherically symmetric) massive body. The interval 5 and the time t 
are measured in time units, r is in length units, c is the local velocity of light, and the constant m is 
defined by 

GM 

m= , ( 2 ) 


1 



where G is the constant of gravitation and M is the mass of the central body. When M is the mass of 
the Sun, 

m= 1.477 X 10 s cm. (3) 


If the central body is the Earth, then 

m = 0.443 cm. 


(4) 


It is well known that in relativity there is no unique system of coordinates in a given space-time. 
The Schwarzschild space-time for a given central body is characterized by one and the same value of 
m, whatever coordinate system may be employed. Two of the simplest alternative systems to the one 
used in Equation 1 are the isotropic and the harmonic systems (Reference 3). In the isotropic system 
the coordinate r is replaced by 7 where 



in terms of which Equation 1 becomes 

, /l - m/2r\ 2 , (1 + m/2F) 4 (dr 2 + 7 2 d£l 2 ) 

ds l = l — dr 

\l + m/2r / c 2 


The analogues of rectangular coordinates may be introduced by 

£ = 7 sin 6 cos 0 , 
tj = 7 sin 0 sin <p , 
f = 7 cos 0 , 

and then the metric becomes 


/l - m!2rV df2 _ (1 +m/2r) 4 (rf£ 2 +rfr? 2 +d$ 2 ) 

\1 + m/2r J c 2 


(7) 


( 8 ) 


The harmonic system has a radial coordinate p where 

r = p + m 


in terms of which the metric has the form 

ds i - 1 ~ m !p dt 2 _ [(1 + m/p)/( 1 - m/p)] dp 2 +p 2 (l +m/p) 2 d£L 2 
1 + m/p c 2 


(9) 

( 10 ) 


Rectangular coordinates somewhat similar to those in Equations 7 are possible with reference to this 
harmonic system (Reference 4). 


To deal with these three radial coordinates at once, it is advantageous to write the metric of the 
Schwarzschild space-time as 


ds 2 = e 2x dt 2 - 


e 2 * dx 2 + e 2v x 2 dti 2 


(ID 
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where X, p, and v are functions of the radial coordinate x involving the constant m. The forms of these 
functions are identified by comparing Equation 1 1 with Equations 1, 6, and 10, respectively, asx is 
identified with r, r, or p. It is clear that, for any one of these identifications, the condition m = 0 
makes X, p, and v zero. Gravitational effects are then absent and the three forms of the metric reduce 
to 


ds 2 = dt 2 - 


dx 2 + x 2 d£l 2 


.*.± 12 * 1 <, 2 ) 
c 2 

with* =r=r = p. The metric is now that of special relativity and the coordinate system (t, x, 6, <f>) has be- 
come the inertial system in which x = 0 is the center of the now massless Sun. This is also an approximate 
description of the coordinate system ( t , x, 6 , 0) employed in Equation 1 1 in those regions of the Schwarz- 
schild space-time in which x is very large compared with m. 

Figure 1 shows portions ABE and FC of the orbits of E and F in space. The Earth is E and the target 
F is another planet or a space probe. It is assumed that the paths are coplanar and lie in the plane 6 = ir/2. 
Figure 1 also shows the paths A MC and CM'B of the outgoing and returning radar signals. They are null 
geodesics of the space-time of Equation 1 1, and it is assumed that the signal is immediately reflected 
at the target on its arrival there. Figure 1 refers to a configuration in which E and F are on opposite 
sides of the Sun. Figure 2 illustrates the configuration when E and F lie on the same side of the Sun. 


Figure 1— Diagram of the paths/WC and CM ' B 
of the radar signals sent between the Earth E 
and the target F. The center of the Sun is at 
O, and OA =x £ , OB = x' E , OC = x p , and OM = X. 



E 



Figure 2-Diagram for the paths AC and CB of the radar signals 
sent between E and F. The paths have been extended to M and 
M', respectively. Other definitions are as in Figure 1. 
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Because E and F are “material particles,” their invariant proper times are nonzero. Consider E 
which is at (x E , ir/2, <p E ) at time t £ and at (x E + dx E , n/2, <t> E + d$ E ) at time t E + dt E . By applying 
to Equation 1 1 the argument that produced Equation 1.8 of Reference 2, it follows that the invariant 
velocity of E is V E , where 

V 2 = [e~ 2 \e 2fl x 2 + e 2v x 2 4> 2 )] E . ( 13 ) 


The dot denotes d/dt and the subscript E means that every term within the bracket must be evaluated 

at x = x E , 0 = 4> e . Also by Equation 1 1 the two events in question are separated by the orbital proper 

time interval . ... 

/ E 2 \ 1/2 

k'7) d, °- (14) 

c ' J E 


ds E = 


As in Reference 2, this terrestrial proper time will be identified with the time kept by atomic clocks on 
Earth. Corresponding formulas for F at C{x p , ir/2, 4> P ) at time t p are 


V p = [e 2 x (e 2ti x 2 +e 2 v x 2 <j> 2 )] p 



F 2 \ 1 / 2 ] 



(15) 


DOPPLER FREQUENCY FORMULAS 

It will be assumed that the radar signal consists of a pulse of electromagnetic waves all of the same 
frequency v. It will be necessary to have a prescription for the frequency relative to clocks that read 
proper times proceeding at different rates. That is, if s and s' are two of these proper times, then 
ds/ds' 1 . The prescription to be employed is that of Kruger. 1 He argues that, in radar experiments, 
it is the phase of the signal which is measured. Let $ be the phase of the wave, and let ds and ds be 
the proper time intervals during which the same increment of phase d$> is measured. Then the fre- 
quencies are to be defined as the rates of change of phase with respect to proper time. Thus, 

v = 

I _ 

v - 

and therefore 

v' 

~ ~ . I • '* iv ' 

v ds 

To express this result in terms of the ( t , x, 0,0) coordinates of the Earth and the target, the 
method described in Section 5.4 of Reference 1 will be employed. This method consists of defining 
the emission and the receipt of a wave by pairs of events as follows. In Figure 1 , let the Earth E be at 


ds ’ 

d7 ’ 

ds 


1 B. Kruger: The Doppler Equation in Range and Range Rate Measurement. NASA GSFC Report X-507-65-385, 1965. 
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A (x p , 7r/2, 4> e ) at time t £ . Then the emission of an electromagnetic wave of frequency v £ will be 
described by the two events E l ( t E , x E , it 12, <fr E ) and E 2 ( t E + dt E , x E + dx E ,tt/2,4> e + d<p E ), where 
dt E may be regarded as the period of the wave in coordinate time. Similarly the receipt of the wave 
of frequency v p by the target Fat Cis defined by the pair of events F 1 (t p , x F , ir/2, <p F ) and 
F 2 (t p + dt p ,x p + dx p , n/2, 4> p + d<p p ). The events (E y , F l ) lie on a null geodesic in the plane 6 = 7r/2 
of the space-time of Equation 1 1, whose differential equations may be written (by Equations 2.308 
and 2.309 of Reference 1) as 



where p is a nonzero parameter varying continuously along the null geodesic and D and )3 are constants 
of integration. These three equations reduce after some calculation to 


d<p _ e 2(x v ^Dc 
dt x 2 


\dt I \ v 2 ) 


It is convenient to introduce here two functions of x: 


QHx,D)= 1 - 


e 2(\- v ) D 2 


(17) 

(18) 


(19) 


P(x,D) = e»~ x Q(x,D). (20) 

The point M on a null geodesic such as AMC of Figure 1 is defined to be the point where x has its 
minimum value, so that dx/dt = 0 at M. Thus if x = X at M, Equation 1 8 shows that 

D = Xe v( - x) - nx) . (21) 


Along the portion A to M of the null geodesic, x decreases as t increases, whereas x and t increase to- 
gether from M to C. Thus the negative square root of Equation 1 8 applies from A to M and the posi- 
tive one from M to C. Hence integration of Equation 1 8 gives for the moments of emission at A and 
receipt at C ( t E and t p , respectively) the relation 


p tg) 



* dx 
Q{x,D) 


( 22 ) 
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Similar arguments apply to the result of dividing Equation 1 7 by the appropriate square root of Equa- 
tion 1 8. Integration then yields 


0/r 0/r D 



e \+y.-2v 

x 2 Q(x,D) 


dx . 


(23) 


Again, the events E 2 and F 2 also lie on a null geodesic whose integrated equations may be obtained 
from Equations 22 and 23 by writing t p dt p 3 t F -t - dtp , Xp dxp , Xp -t dXp , (j)p d(pp , and 
(pp + d(p F for t E , tp, x £ , x F , <P E , and <p F , respectively, and also by altering the constant of integration 
D to D + dD. It is also easy to show that 

/ D dD e 2(x-y) \ 

Q~Hx,D + dD) = Q~Hx,D) (1 + . (24) 

\ x 2 Q 2 {x,D) ) 


The differences dt p - dtp and d<p F - d<t> E are now required; and, since the second difference probably 
involves the more complicated calculation, it is worked out in detail. 

The null geodesic equation (Equation 23) for E y and F, and its analogue for E 2 and F 2 yield 


d<p P - d<p F - (D + dD) 


Xp+dXp 





p \+H-2v 

— ! dx 

x 2 Q(x, D + dD) 


-D 



+ 



e \+ P -2v 

x 2 Q{x,D) 


dx , 


where dX is obtained in terms of dD by taking differentials of both sides of Equation 21 . If the right- 
hand side of the last equation is evaluated to the first order in the differentials, it is found after some 
calculation and the use of Equation 24 that 


D 

( d<t> F — d<Pp ) — 
c r c 


“(f 



e \+H-2v 

x 2 Q 3 (x,D) 

D 2 

+ — [J{x E ) dxp + J{x F ) dXp 


2J(X) dX) , 


where, in the last three terms, 


A corresponding calculation applied 
D dD i 


dt F dt p 


r* E r* 

\Jx Jx 


p \+V-2v 

J(X) = . 

x 2 Q(x,D) 

to Equation 22 and its analogue for F 2 and F 2 yields 

p \+n- 2v 

dx 

x 2 Q 3 {x,D) 

+ - [I(x E )dx E + I(x F ) dxp - 2 1(X) dX] , 


(25) 


(26) 
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where 


/(*) = 


e n-\ 

Q(x,D ) ' 


If now Equation 25 is subtracted from 26 and Equations 19 to 21 are also used, it follows that 


D dXp D “•* e 

dt F - — - P(x F , D) = dt E d(j> E + P(x £ , £>) 


dx E 


( 27 ) 


which is the final result for the outgoing signal A to C. It has been established without the necessity of 
finding dD or dX. Presumably Equation 27 is still valid when x f = X so that C and M coincide, in spite 
of the fact that in these circumstances Q(X, D) = 0 and therefore also P(X, D ) = 0. The result is also 
expressible in terms of the Mime rates of change ( dx/dt , d<j>/dt = x, <j>) of the coordinates of E and F. 

In fact, Equation 27 is also 


dt E 1 - D4> F /c - P(x F ,D)x F /c 
dt F 1 - D<i> E /c + P(x E ,D)Xp/c 


(28) 


The returning signal is immediately reflected at C so that none of the symbols with subscript F 
are altered. But the point of closest approach is now M' in Figure 1 , for which the constant is D' . The 
Mime period is dt' E at B whose coordinates and their rates of change are x ' E , <p ' E , x ' E , and <p' E . The sig- 
nal is also now traveling from C to B. When all these changes are taken into account, and the preceding 
calculations are repeated, it turns out that 

dtp 1 +D’4>' E /c- P(x' D')x'/c 

= t (29) 

dt' E 1 +D'<pp/c+P(Xp,D')Xp/c ' 


Formulas 28 and 29 refer to the configuration of Figure 1 . When C, A , and B all lie “on the same 
side” of the Sun, as in Figure 2, the calculations are simplified because, in the integrals corresponding 
to Equations 22 and 23, the Emit X does not occur. The limits in the integrals are simply x E and x F . 
The pathv4C of the outgoing signal now yields 


dt E 1 - D(j>p/c + P(Xp, D)x f /c 
dtp 1 - Dipp/c + P(x E , D)Xp/c 

while, for the. returning signal along CB, 


dtp 1 + D'<j)'p/c - P(x'p , D’)x' e /c 
dtp 1 + D'<j>p/c - P(Xp,D')Xp/c 


(30) 


(31) 


The final step consists in transforming from coordinate-time intervals to proper time intervals and 
thence to frequencies through Equation 1 6 and results such as Equations 28 to 3 1. A frequency v E 
emitted from A is received as v p at C, reemitted as v p at C and received as v' E at B. The two-way radar 
Doppler frequency formula is then 

v'p ds E [e*(l -V 2 /c 2 ) 1 I 2 ) e dtp 

— = = , (32) 

v E ds'p [<? x (l - V 2 /c 2 ) 1 I 2 ] e , dt'p 


1 



and the one-way radar formulas are 


V F d$E 

[e x {\ - V 2 /c 2 ) ll2 ] E dt E 

V E 

ds F 

[e x (l - V 2 /c 2 ) l ' 2 ) F dt F 

v 'e 

ds F 

[c x (l - V 2 /c 2 ) 1/2 ]p dtp 


v F ds’ E [e*(l- V 2 /c 2 ) l ' 2 ] E , dt' E ‘ 


( 33 ) 


(34) 


In these formulas the ratios of the coordinate-time differentials are obtained from Equations 28 and 29 
for the Figure 1 configuration and from Equations 30 and 31 for that of Figure 2. 


Though Formulas 33 and 34 have been called Doppler frequency formulas, they nevertheless 
incorporate the gravitational red-shift effect. Consider the case for which E and F are at rest, so that 
in Figure 1 , for example, A and B are coincident and x' E = x g . Though now x, <j>, and V are all zero at 
C and A, yet v p fv E and v' E /v F are in general different from unity because \(x F ) =£ M* E )- However, 
in the two-way radar formula (Equation 32) the gravitational red-shift effect has been eliminated and 

v 'e =v e- 


SPECIAL RELATIVITY 

The results of the preceding section may be elucidated by considering their meanings in the 
degenerate case of special relativity. In this theory, X, n, and v are all zero, the geometry is Euclidean, 
and the paths of the signals are Euclidean straight lines. Hence by Equation 1 3 the velocity of either 
E or F is given by 

V 2 = x 2 +x 2 0 2 

= u 2 +v 2 , (35) 

and, by Equation 21 , D = X and D' = X'. The angle \p may be defined by 

— = cos y}j , (36) 


and then Equations 1 9 and 20 show that 

/ D 2 \ 112 

Pi x, D ) = Q(x, D)= ( 1 —J = sin $ . (37) 

Replacement of D by D' in the last two formulas gives corresponding results for an angle \p'. The 
velocity diagram for the outgoing signal, A to C, of Figure 1 is shown in Figure 3. The one-way 
Doppler frequency formula (Equation 33) may be converted by Equations 28 and 35 to 

1/2 

1 - V 2 lc 2 \ 1 - ( D/x f )(v f /c ) - P{x f ,D)u f /c 


(1 - V 2 /c 2 1 - (D/x f )(v f /c) + P{x F , D)u f /c 
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Figure 3— Velocity components of E at A and F at C. 


But 


— v f + P(x f ,D)u f = u F cos \p F + u F sin \jj F 
x F 


— v £ - P(x e ,D)u e = v e cos \p E - u E sin \jj E 
X E 


where the velocity components w p and w E , of F and E along the (straight) line of sight, are reckoned 
to be positive in the direction A to C. Thus 


V F (l - v \l c2 ) 112 1 ~ w F /c 

V E (1 ~ Vj,/c 2 ) 112 1 - w E /c 


(38) 


If F is identified as an “observer” who is at rest at C so that w f = 0 and V p = 0, then the emission 
frequency in the rest frame of E at A is v £ and the frequency received by the observer at C is v F . In 
terms of wavelengths, X £ - c/v E and X F = c/v F \ thus the last equation is 

X F 1 + ( ~w E /c ) 

— = ~ , (3 

X E (1 - V E /c 2 ) 1 ^ 2 


where the case of a negative value of w £ is presupposed so that E has a velocity of recession along *he 
line of sight as viewed from C. Equation 39 is the well-known Doppler formula of special relativity 
(Reference 5) including the transverse Doppler effect represented by the term in V E /c 2 . 

The velocity diagram for the returning signal C to B is shown in Figure 4. The factors linear in 
the velocity components are those shown in Equation 29. The w components are again taken in the 
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Figure 4— Velocity components of E at B and F at C. 


direction of motion of the signal. Hence 

D' 


v F +P(x F ,D')u F = v p cos \p’ F + Up sin \jj' F 


= - w. 


D' 

— v' F - P(x'D')u' = v' cos i //' - u' sin \jj' E 


= - w 


E ■ 


Thus the special relativity version of Equation 34 becomes 

^1 - V 2 /c 2 \ l/2 1 - w'/c 


a - v' 2 !c 2 l - w F /c 


(40) 


which is essentially the same formula as Equation 38 when it is remembered that v p now plays the 
part of the emitted frequency. The reader may draw his own velocity diagrams for the Figure 2 con- 
figuration. It is to be borne in mind that the v components are drawn in the direction of <t> increasing 
and the w components are always in the direction of motion of the signal. Then Equation 38 again 
follows from the special relativity versions of Equations 30 and 33, and Equation 39 follows from 
those of Equations 31 and 34. 

The main point of this discussion is that the ratios of the Mime differentials in Equations 32 to 
34 give the effects of the line-of-sight velocities of E at A and B and F at C and this interpretation may 
be expected to hold for the non-Euclidean geometry of the Schwarzschild space-time. In this case 
the principal complication arises through the definitions of the angles \p and i//. An illustration, based 
on the use of the metric of Equation 1, will be found in Section 3 of Reference 2. 
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NUMERICAL VALUES 


To obtain an estimate of the orders of magnitude involved in the terms of the general relativity 
formulas^ the target will be assumed to move in the solar system under the gravitational attraction of 
the Sun alone, which is in accordance with the definition of the Schwarzschild space-time. Inspection 
of formulas 28 to 34 shows that the arithmetic will be simplified if the Earth is regarded as at rest 
during the radar experiment. Under these conditions, the Figure 1 configuration has A and B coin- 
cident, so that D = D'. Then Equation 32, with the aid of Equations 28 and 29, reduces to 


v' 1 - Dfi/c - P(x, D)x/c 

— = ; (< 

v £ 1 + D<t>/c + P(x, D)x/c 

where the subscript F has been omitted on the right-hand side. The expression to be evaluated may 
then be written 


V(x) = 


D (x4>) 
x c 


+ P(x, D)- . 
c 


(42) 


Suppose that it is decided to use the coordinate r of Equation 1 in the first instance. 


Equation 1 1 and 


e 2x _ 


2m 



Hence x = r in 


so that Equation 2 1 with X — R is 

and by Equations 19 and 20 

P(r, D ) = 



e 2v = 1 , 


D 




2mV 1/2 

1 “ T/ 


2m 

r 




j R 2 (1 - 2 m/r) 
r 2 (1 - 2 m/R) 


1/2 


(43) 


(44) 


(45) 


These are exact formulas that would have to be used if r were approximately equal to R ; i.e., for a 
configuration in which C lay close to M. In a large number of configurations, however, the ratio 
R 2 /r 2 is substantially less than unity. Thus it is possible to use expressions for D and P that are of the 
first order in m/r and m/R, namely, 


D = R 



P(r, D) = 


' * 2 V ' 2 
: ■ 7 2 ) 


2m/ r- (m/r + m/R)(R 2 /r 2 ) 
1 - R 2 /r 2 


(46) 
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The exact expression for r\{r) is obtained by substituting x = r and using Equations 44 and 45 in 
Equation 42. The approximate expression derived from Equation 46 is 


V(r) = 


Mr) 

c 


1 + 


m r<p 
r Mr) 


2 m/r - (m/r + m/R)(R 2 /r 2 ) r 


(1 - R 2 /r 2 ) 112 


w(r) 


where 


R 


Mr) = ~ ( r<t> ) + 

r 




( 47 ) 


(48) 


In Formula 47 and in all subsequent equations for 17, the contributions of the terms involving m (the 
m terms) are enclosed in square brackets with subscript m. It was shown that when m = 0, and special 
relativity applies, r? reduces to w, which is the radial velocity of F at C relative to the stationary Earth at A 

General relativity, however, contains no prescription that compels the use of r in preference to F, 
to p, or to any other radial coordinate that may legitimately be used in the Schwarzschild space-time. 
Suppose, therefore, that x is next identified with F so that now 

\1 + m/2r/ 


An alteration in the value of m is, of course, forbidden because the same central body is in question 
Then Equation 2 1 with X = R gives 

- -(1 + m/2R) 3 


D = R 


and Equations 1 9 and 20 yield 

- (l+m/ 2 F ) 3 

P(r,D)= 

1 - m!2r 


1 - 


1 - m/2R 

£ 2 (1 +m/2R) 6 (1 - m!2r) 2 


r 2 (l - m/2R) 2 (\ + m/2r) 6 


1/2 


(49) 


(50) 


The approximate form of 17(F) derived from these expressions and from Equation 42 is then 


77(F) = 


Mr) 


1 + 


2m r<p + 2 [m/r- (m/R)(R 2 /r 2 )] r 


w(F) 


(1 -R 2 /r 2 ) 112 w(F) 


(51) 


where Mr) is obtained from Equation 48 by replacing R and rby R and r, respectively. Lastly the 
coordinate p of Equation 10 could be employed so that x = p and 


e 


2 \ 


£ 2 M 


e 


2v 


1 - m/p 
1 + m/p 

1 + m/p 
1 - m/p 
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Then, if X = <R and D = A, Equations 19 to 21 yield 


(1 + m/<R) 3/2 

<K 

(1 - m/<R) 1/2 


Pip, A) = 


1 + m/p j (R 2 
1 - m/p p 2 


(1 + m/<R) 3 

1 - m /( R 


1 - m/p 
( 1 + m/p) 



(52) 

(53) 


It is noticed, by comparing Equation 49 with 52 and Equation 50^ with 53, that the exact dependence 
of D on R is different from that of A on (R and similarly for P(r, D ) as compared with P(p, A). But to 
the first order for the m terms 



2m 2 [m/p - (m/(R)((R 2 /p 2 )] p 

P vv(p) (1 - <R 2 /p2)i/2 W (p) 


(54) 


where w(p) is obtained from Equation 48 by writing (R and p for R and r, respectively. Comparison of 
Equations 51 and 54 shows that r\{r) depends on r and R in the same way as r?(p) does on p and <R. To 
this order of approximation therefore, it is immaterial whether the isotropic coordinate r or the 
harmonic coordinate p is employed. But this would not be true for configurations in which r were 
nearly equal to R and p to <R. 

The Newtonian analogue for the motion of F is given by the elementary two-body problem in 
which the mass of the Sun M e is regarded as infinitely great compared with the mass of any planet or 
space probe. As an example, let F be a space probe whose orbit, calculated by Newtonian theory, is an 
ellipse with the Sun at one focus, of perihelion distance^/ 5 and eccentricity e n =2/3, where a n is the astro- 
nomical unit of Newtonian orbit theory. The semimajor axis of the orbit of F is then 2a J 5 and the semi- 
latus rectum is a n /3. At any instant of absolute Newtonian time t n , the radial and cross-radial components 
of the velocity of F are 


dr„ 


dt„ 


d<P„ 


v = r — 
" " dt 


where (r ,<f> ) are the coordinates of F computed by Newtonian gravitational theory. They are 
coordinates in Euclidean space; such Euclidean coordinates, which arise from Newtonian theory, will 
be called N coordinates. Now let Figure 1 refer to a configuration in which F is at that end of the 
semilatus rectum of its orbit at which r n (equal to a n /3) and <j> n are increasing together. The choice 
is dictated by the well-known Newtonian theorem which states that the radial component of velocity 
is greatest at an end of the latus rectum. At such a point 


n 

C 


(GM P 


0 



(55a) 
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1/2 


(55b) 


IGM, 


o 


1 - e n 3a n, 


Incidentally, reference to Equation 2 shows that both of these are of order (m/a) 1 ^ 2 , a consideration 
that might be employed in what follows. However, it is probably easier to go directly to the numer- 
ical values, which are given in Equations 56. 


In addition to placing F at one end of the latus rectum of its orbit, it will be assumed that F is 
seen from E at grazing incidence to the solar surface. In the sequel this will result in the m terms of 
Equations 47, 51, and 54 being as large as possible. The r n coordinate of the point M of closest 
approach to the Sun’s center will be identified with the solar radius R Q . 

It is assumed that the quantities in the right-hand sides of Equations 47, 51, and 54 may be 
identified with the corresponding Newtonian values. Thus there are three alternative identifications 
according to whether * is chosen to be r, r, or p, and X to be R, R, or <R. The three alternatives can be 
represented by (Reference 6). 



= ^(1.496 X 10 13 ) cm, 
X = R q 

= 6.960 X 10 10 cm, 

1 = U _n 

c c 

= 1.147 X 10" 4 , 

X A = V JL 

c c 

= 1.720 X 10“ 4 , 
gm q 

m = 

c 2 

= 1.477 X 10 5 cm. 


(56) 


The values of u and v 

n n 


are calculated from Equations 55 with e n = 2/3. It then follows from 
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Equation 48 and its analogues for the coordinates r and p that 

w{r ) _ w(r ) 
c c 


w(p) 


c 


and that the m 
respectively, 


- 1.171 X 10 - 4 [ 1 - 9.539 X 1(T 5 +0(1(T 9 )] , 
terms, if r? is to be correct to the order 10 “ 12 in Equations 47, 5 1, and 54, are, 

[r] m = 1.011 X 1(T 7 , 

r 

= 1.443 X 10 - 7 . 


( 57 ) 


(58) 


Thus, for the particular configuration in which F is at one end of the latus rectum of its orbit and is 
seen from E at grazing incidence to the solar surface, the m terms differ in value by as much as 42 per- 
cent as r, r, or p is identified with the N coordinate r n . This is not a property peculiar to this con- 
figuration. Examination of the m terms in Equations 47, 5 1 , and 54 shows that differences would be 
found in other configurations also. 

A two-way radar experiment would, of course, give v' E /v E rather than 77 . However, Equation 41 
to a sufficient accuracy is 

v j_ - 

V E 1 +r * 

= 1 - 2 r 7 + 2r) 2 . 

Thus for the r coordinate, Equations 57 and 58 give 
v 'e 

— = 1 - 2.342 X 1CT 4 +4.978 X 1(T 8 + 0(1(T 12 )- [2.368 X 10 " 11 +•••],„ - (59) 

V E 

If either r or p had been employed instead, the amount of the leading m 
[3.379 X 1 0“ 1 1 l m • This example shows that an accuracy of one part in 
measurement of v' E /v E if the values of the m terms were to be deducible 
experiments. 

The motion of the Earth produces additional mathematical complications in the evaluation of 
v' e !v e from Equation 32 and either Equations 28 and 29 or 30 and 31. One complication arises be- 
cause the difference between D and D' must now be calculated. Another is that the values of these 
constants are not known directly, as has been assumed in the foregoing example. In fact only the 
coordinates of A, C, and B can be regarded as calculable from the orbits of E and F. The quantities 
D and D' must therefore be expressed in terms of the coordinates of these three points. An approxi- 
mate method of dealing with these matters, in terms of the r coordinate, is given in Reference 2. 


term would have been 
1 0 1 1 would be required in the 
from two-way radar 
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COORDINATE IDENTIFICATIONS 


The smallness of the m terms and the considerations mentioned at the end of the last section 
indicate that orbits of high accuracy for E and F are needed so that the coordinates of A, C, and B 
and the velocity components at those points can be predicted as exactly as possible. To obtain such 
orbits, it is necessary to allow for the attraction of the planets on each other or on a space probe. This 
is done for the American Ephemeris by a method which is “in accordance with the Newtonian law of 
gravitation, modified by the theory of general relativity” (Reference 7). This means that the Euclidean 
space and absolute time of Newtonian theory are employed in writing down the equations of motion 
of planets and space probes. These equations contain terms that allow for the motion of the Sun and 
for the mutual attractions of the planets. But in addition, the equations contain terms, essentially 
ad hoc, that produce a small rotation of the line of apsides of a planetary orbit. This small rotation 
has the same value as that predicted for the orbit of a test particle in the Schwarzschild space-time, 
where, as already stated, the planetary perturbations are by definition zero. A similar procedure 
appears to be employed at the Lincoln Laboratory (Reference 4). A relativity perturbation is calcu- 
lated for a test particle in the Schwarzschild space-time in terms of harmonic coordinates. It is 
multiplied by an arbitrary constant factor (to be found eventually from observation) and is then 
inserted into the complete Newtonian equations of motion to produce an additional perturbing force. 

In either case, it can be argued that N coordinates, and their associated Euclidean velocity components, 
are obtained. Both methods will be called “classical celestial mechanics.” 

The problem of identifying the N coordinates with some set of coordinates in the non-Euclidean 
space of the Schwarzschild field has still to be solved. One way of proceeding is by definition, as was 
done explicitly in the preceding section. Or the definition may be made implicitly as in the Lincoln 
Laboratory method. The insertion of the relativity perturbation into the Newtonian equations of 
motion carries with it the implications that the t time of the Schwarzschild field is identical with the 
Newtonian absolute time t and that the harmonic coordinates are identical with the N coordinates. 

It is true that Fock (Reference 3) has attempted to prove theoretically that harmonic coordinates are 
in some sense superior to others, a contention that contradicts Einstein’s aim of putting all coordinate 
systems on a par with one another. Fock’s idea might carry weight if it could be shown that the m 
terms in the rj(x) of Equation 42 always had their smallest numerical values when harmonic coordinates 
were employed. It might then be argued that these coordinates most closely approximated the N 
coordinates. But it has been shown that harmonic coordinates do not always produce the smallest 
values of the m terms. 

When radar time delays were considered (Reference 2), the m terms were found to be larger, of 
the order of 1 0~ 7 , than those shown in the Doppler frequency formula (Equation 59). They are also 
different in value according to the identification employed. Thus the same problem of coordinate 
identification arises. 

An alternative solution of the problem is possible in principle. Orbits of E and F would be 
calculated by classical celestial mechanics combined with observations of E and F for those configura- 
tions in which the m terms could be expected to be negligibly small. The N coordinates thus produced 
would then be identified in turn with the coordinates used in Equations 1, 6, and 10 to predict the 
values of the m terms for those configurations in which these terms were expected to be as large as 
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possible. Simultaneously, measurements would be made of time delays and Doppler frequency ratios 
in these configurations with large m terms. In this way, not only would the predictions of general 
relativity be checked, but the Schwarzschild coordinate system that most closely corresponds to the 
N coordinates would be identified. 

In conclusion, reference may be made to a point that is more easily discussed with respect .to 
Doppler frequency formulas than to time delays. The frequency ratios, in Equation 16 or 32, for 
example, are obtained as the ratios of two invariant intervals ds. It might therefore be argued that the 
particular coordinate system used in the calculations could not affect the result. This is true so long 
as coordinate transformations such as those in Equations 5 and 9 are in question. In these transforma- 
tions, to every r there corresponds a different r or p calculated from the two formulas, and vice versa. 
But this is not the operation envisaged in the coordinate-identification procedure. In the previous dis- 
cussion, for example, the question was whether r, r, or p should be identified with a n l 3; this is not 
a transformation of the kind envisaged in Equations 5 and 9. Thus, invariance of the Doppler 
frequency ratio is not to be expected. 
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